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ABSTRACT 

In this work we present a reactions module for "Chaos Many-Body Engine" (Grossu et al., 2010 [1]). 
Following our goal of creating a customizable, object oriented code library, the list of all possible reactions, 
including the corresponding properties (particle types, probability, cross-section, particles lifetime etc.), could 
be supplied as parameter, using a specific XML input file. Inspired by the Poincare section, we propose also the 
"Clusterization map", as a new intuitive analysis method of many-body systems. For exemplification, we 
implemented a numerical toy-model for nuclear relativistic collisions at 4.5 A GeV/c (the SKM200 
collaboration). An encouraging agreement with experimental data was obtained for momentum, energy, 
rapidity, and angular ri distributions. 
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1 



Restrictions: 



Additional comments: Easy copy/paste based deployment method. 
Running time: quadratic complexity 



1. Introduction 

In a previous work [1], we presented "Chaos Many-Body Engine" code library for chaos 
analysis of relativistic many-body systems. Inspired by existing studies on Fermi nuclear 
systems [2-5], we tried to apply chaos theory to the more complex case of nuclear relativistic 
collisions at 4.5 A GeV/c. In this domain of energies it is possible to apply the many-body 
representation, as the nucleon de Broglie wavelength is less than the average inter-nucleonic 
distance, and the nucleon mean free path is shorter than the target radius [6-10]. On the other 
hand, the two-nucleon reactions are a key point in understanding the processes involved in 
nuclear collisions. In this context we motivate our efforts of developing a C# [11] reactions 
module, integrated with "Chaos Many-Body Engine". 

2. Two particles reactions and decays 

The second version of "Chaos Many-Body Engine" library is designed as a general 
numerical solution for the simulation of three-dimensional relativistic many-body systems 
with reactions. We implemented: two particles reactions, decays, and stimulated decays. 

a+b->c+d reactions. For simplicity, we considered that only the momentum component 
oriented towards the direction between the two bodies is changing during the collision. 
Working in a Cartesian system in which this direction corresponds to abscise, by applying the 
relativistic energy and momentum conservation laws, wrote in the natural system of unities, 
we obtained the following second degree equation: 

(a 2 - 1)E 2 + 2 « /?£ c + /? 2 + p 2 a + p 2 a + m 2 c = (1) 

with: 

E a+ E b 

a = (2) 

Pxa "I" Pxb 

™ 2 a + Py b + pit - {E a + E b ) 2 + (p xa + p xb ) 2 - p 2 a - p 2 zb - m 2 

2(Pxa+Pxb) 

where, m a , m h , m c , m d are the rest masses, p xa , p x b, Pxc, Pxd, p ya , p y b, p yc , p y d, p za , Pzb, p zc , Pzd 
the Cartesian components of momentum, and E a , Eb, E c , Ed are the particle energies. 
The momentums of the two resultant particles are given by: 

Vc = ±^E? ~ m 2 c (4) 
Vd=Va+Vb~ Vc (5) 

Only the physical solutions of equations (1-5) are chosen. We impose also two additional 
conditions: 

The two colliding particles must approach each other 
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The distance between the two particles must be less than a specific value, r max , 
calculated in agreement with the cross section a: 

a->c+d decays. For simplicity, we assumed that each decay take place after the 
corresponding lifetime, measured in particle's own reference system. Thus, the program was 
improved in order to adjust each particle lifetime, according to the corresponding 
instantaneous velocity: 

dt 

t 2 = t 1 + ; (7) 
where t is the time, /?=v/c, and dt is the integration interval. 

In the frame of decaying particle, the two new particles are emitted on a random 
direction, with the following momentum: 

1 , 

V =Pc=Pd = ~ VN - (rn c + m d ) 2 ] [m a - (m c - m d ) 2 ] (8) 

2m a 

For obtaining the result in the simulation reference frame, we must relativistic compose 
the two calculated velocities, with the transport velocity. 

Stimulated decays are decays triggered by collisions. At least one participant must have 
at least one decay schema. This kind of reactions is implemented by simply setting the 
particles lifetime values to zero on collision events. 

It is important to notice that more complex schemas could be considered by employing 
various combinations of the, previous described, basic reactions. For example, for a + b -> c 
+ d + e, one can use a virtual particle "f ', with zero lifetime, which decays into "d" and "e". 
Thus, the engine will process the following sequence: a + b->c + f;f->d + e. 



3. Program description 



Our main goal was to create a highly configurable reactions module, integrated with 
"Chaos Many Body Engine". Thus, the list of all reaction schemas which applies to the 
many-body system of interest is stored in a DsReactions typed dataset object. The main 
advantages of using datasets [12] are related to XML serialization, and very easy filtering 
capabilities, based on SQL-like query strings. 

The DsReactions dataset contains three datatables. The DsReactions. Particle table is used 
for defining all necessary particle types. It contains the unique identifier (IdParticle), and the 
most common particle properties (mass, charge, and lifetime). The 
DsReactions. ParticleProperties datatable is used for dynamically defining particle properties 
(spin, isospin, color charge and so on). It is in a many to one relation with 
DsReactions. Particle, and contains each property name and value. One new Particle class 
instance can be created, based on information stored in the DsReactions. Particle and 
DsReactions. ParticleProperties datatables, using the NBody.NewP article method, which 
expects the particle type {IdParticle) as parameter. 
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The DsReactions. Reactions datatable is used for defining the list of all possible reactions 
for the specific system of interest. The IdParticlelnOl and IdParticleIn02 columns represent 
the unique identifiers of the input particles, while IdParticleOutOl and IdParticleOut02 are 
the reaction output identifiers. The MaxDistance column is directly related to the cross- 
section parameter (Eq. (6)). For avoiding also any unwanted effect related to infinity values 
near origin, we considered a minimum distance (the MinDistance column), under which the 
reaction is forbidden. For decays, this property has a different signification. It represents the 
initial distance between the two, new generated, particles. When more channels are possible, 
the reaction probability can also be specified. The same datatable is used for all reaction 
types. Thus, a decay is specified by setting IdParticleIn02=null, while for a stimulated decay, 
one will set both IdParticleOutOl , and IdParticleOutOl to null. 

The list of all allowed reactions can be supplied/modified by simply creating/changing the 
XML file used as data source for the DsReactions dataset. The user should carefully define 
each reaction, in agreement with all specific conservation laws (charge, baryonic number and 
so on), and avoiding also circular references. 

Based on a second order Runge-Kutta algorithm, the Nbody.Next method is used for 
"moving" the many-body system into its next status [1]. We improved this method for 
finding also all particles susceptible to participate to reactions. Thus, we fill one generic list 
of Reaction objects with all constituent pairs for which the mutual distance is less than the 
maximum cross-section radius, Eq. (6), and another one, with all particles for which the own 
time, measured in the own reference system, is greater or equal with the corresponding 
lifetime. The two lists are used as input for the reactions engine. For optimization reasons, 
the reactions module is called with a frequency which can be specified as parameter. It is 
important to notice that the program will automatically increase this frequency in certain 
circumstances (the number of particles susceptible to participate to reactions changed, and/or, 
at least one reaction took place in the current iteration). 

The ReactionEngine class contains the main processing routines. The GetReactions 
method assures that one constituent will participate to only one reaction. The decays are 
treated in priority. The CheckReaction method is used for finding all reactions which apply to 
the input data (the, previous described, generic lists), according to the information stored in 
the DsReactions dataset. The GetReactionOutput and GetDecayOutput methods select only 
the physical solutions which satisfy the energy and momentum conservation laws (equations 
(1-5)). If there are more possible reaction schemas for the same input, one solution is 
randomly chosen, according to its probability (the DsReactions. Reactions. Probability 
column). 

The Nbody.ProcessReactions function is used for applying the processing result to the 
NBody instance. It removes the old particles, and adds the new ones from/to the mParticles 
generic list. The OnReactions event is raised if at least one reaction took place. As the 
collisions are inelastic, the kinetic and potential energies are also adjusted after each reaction 
event, in order to avoid any impact on the energy conservation precision test [13]: 



-log 
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E(t) - E(t = 0) 



E(t = 0) 



(9) 



The structure of the simulation output files was changed according to the new necessities, 
trying to assure also the compatibility with the previous version. Thus, in the header file we 
added the creation time of each particle (the Time column), an integer (Particle) which 
uniquely identifies each particle instance, and the particle type (Id), which corresponds to the 
IdParticle column defined in the DsReactions dataset. As the composition of the system 
could change in time, the unique instance identifier (Particle) was added also to the data file. 
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5. Chaos analysis of many-body systems 



In [1] we implemented the multi-dimensional Lyapunov Exponent [14,15]: 



1 d(t) 

I«Hm-Jn-^=limI(t) (10) 
t->oo t a(UJ t^oo 

where d(t) represents the phases space distance between the two systems: 



d = 



it, 

I 



(0W; 2 ) 2 + (p £l -p £2 ) 2 ) (11) 



where r tl , p ilr are the position, respectively the momentum, of the particle "f from the 
first system, and r t 2 , pn, are the analogous coordinates of the second system. 

As the structure of one many-body system with reactions could change in time, the 
previous distance cannot be applied any more. We propose instead a "global" approach: 



d = 



N 



+ (l2) 



1=1 



t=l 



The BaseSystemData. GlobalLyapunovExponent static method is implementing the 
Lyapunov Exponent based on this distance definition. 

Inspired by the Poincare section [15,16], we propose also a new intuitive investigation method. 
Thus, for one many-body constituent we considered the following set: 



Mi = {(xtQtUtW) |-|f(t) = V ^i(t) = V-^(t) = o} (13) 
We defined the "Clusterization Map" as: 



M 



= []M i (14) 



i=l 



Using this technique, one obtains a "bi-dimensional projection" of all clusters created 
inside the system. For a better precision, the BaseSystemData. ClusterizationMap static 
method is searching all points for which the derivatives of momentum components signs are 
changing. 
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Fig.l. Two Clusterization Map examples for Cu+Cu at 4.5 A GeV/c nuclear collision simulation, obtained with 
Chaos Many-Body. (Left) Central collision. (Right) Non-central collision. 



6. Application to the C-C nuclear collision at 4.5 A GeV/c 



As an example of use, we implemented a simplified model for nuclear collisions at 4.5 A 
GeV/c (the SKM200 collaboration [8,17-19]). The colliding nuclei are represented as two 
distinct sets of nucleons, initially placed in the vertices of a cubical network, with the edge 

1/3 

calculated in agreement with the nucleus radius r=roA . The target is initially at rest, while 
the incident momentum of the projectile constituents could be specified as parameter. For 
simplicity, the Lorenz contraction is ignored. 

We employed a finite depth Yukawa potential well, together with a coulombian, and a 
short distance repulsive term: 



JJi 

e a 



—k ln(r) , r t j < 1.134Fm 

mi 
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Aner, 



— k ln(r), 1.134Fm < r £; - < 2.2Fm 



(9) 



V 4- qiqj 



Aneri 



rij > 2.2Fm 



where Vo=35MeV is the depth and a=2Fm is the radius of the potential well, ry 
represents the distance between the two bodies, and q is the electric charge. We empirically 
chose k=200. 

The following interactions were considered: 
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p + n -» p + A 
p + n -> A + + n 
p + p -> A + + p 
n + n -> A + n 

A°^ p +7T~ 

A°^n + 7r° 

A + ^ p +7T° 

A + -> n + n + 
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Working with a 10" Fm/c value for temporal resolution, we simulated 500 C+C events, 
with random collision parameters. The comparison with experimental data illustrates an 
encouraging result for momentum, energy, rapidity, and angular n distributions. 




-■— many-body engine 

* exp SKM200 
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Fig.2. Comparison between Chaos Many Body Engine and experimental n distributions for C+C at 4.5 A 
GeV/c. (Up Left) momentum distribution. (Up Right) energy distribution. (Down Left) rapidity distribution. 

(Down Right) angular distribution. 



7. Conclusion 



Starting from exiting analysis of Fermi nuclear systems, one of our main goals was to 
apply chaos theory to relativistic nuclear collisions at 4.5 A GeV/c (the SKM200 
collaboration). In this context, we developed a C# reactions module for "Chaos Many-Body 
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Engine". An important attention was paid to the application flexibility and extensibility. 
Thus, we improved the code in order to accept any number of, dynamically loaded, particle 
properties (spin, isospin, color charge and so on). The list of all particles and possible 
reactions could be also supplied as a simulation parameter using a specific XML input file. 

As the structure of one many-body system with reactions could change in time, in the 
implementation of the multi-dimensional Lyapunov Exponent we considered a particular 
definition for the distance between two systems, equation (12). We propose also the 
"Clusterization Map", equation (14), as a new analysis method, which can bring intuitive 
information on the clusters created over the time in the system. 

The, previous described, example of use (C+C collision) does not require significant 
resources to run (CPU 1.0GHz, 128M free RAM, .Net Framework 2.0 running on MS 
Windows XP or later). However, the resources could become a critical problem with the 
increasing of reactions and particles numbers (e.g. Cu+Cu or Au+Au collisions). 

Using a simplified model for C+C nuclear collision at 4.5 A GeV/c, we obtained some 
encouraging results for momentum, energy, rapidity, and angular n distributions (Fig.2). 
Further analysis along those lines is currently in progress. 
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